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Abstract 

Preliminary results are presented which suggest that scaling and singu- 
larity characteristics of solar wind and ground based magnetic fluctuations 
appear to be a significant component in the solar wind - magnetosphere in- 
teraction processes. Of key importance is the intermittence of the "magnetic 
turbulence" as seen in ground based and solar wind magnetic data. The meth- 
ods used in this paper (estimation of flatness and multifractal spectra) are 
commonly used in the studies of fluid or MHD turbulence. The results show 
that single observatory characteristics of magnetic fluctuations are different 
from those of the multi-observatory AE-index. In both data sets, however, 
the influence of the solar wind fluctuations is recognizable. The correlation 
between the scaling/singularity features of solar wind magnetic fluctuations 
and the corresponding geomagnetic response is demonstrated in a number of 
cases. The results are also discussed in terms of patchy reconnection pro- 
cesses in magnetopause and forced or/and self-organized criticality (F/SOC) 
of internal magnetosphere dynamics. 



1 Introduction 

According to recent knowledge magnetospheric plasma fluctuations may be ex- 
plained in terms of multiscale intermittent models including nonlinear self-organization 
processes (forced or not) near criticality (F/SOC) |Chang(1992)| |Chang( 1999)1 



Chap man et al.(1998)| |Takalo et al.(1999)| |Klimas et al.(2000)| . The basic phys- 



ical concepts and the corresponding path-integral and renormalization group for- 
malisms of nonlinear self-organization of magnetospheric processes are described 
in [Chang(1999) ]. These models predict multifractality, non-Gaussian, power-law 
distributions for certain measurable physical quantities and associated power-law 
scalings for power spectra. In fact, several experimental studies seem to support the 
theoretical predictions quite well [ Consolini et al.(1996)| |Milovanov et al.( 1 996) , 
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Consolini(1997)]|Borovsky et al.(1997)||Consolini and De Michelis(1998)||Voros (1998)1 



Uritsky and Pudovkin (199T)]|Consolini and Lui(1999)HChapman et al.(1999)HFreeman et al.(2000a)i 



Freeman et al.(2000b)| jSharma et al.(2001)| . Scale-free power spectra of magne- 
tospheric fluctuations, however, admit several possible explanations other than 
F/SOC. [Wat kins et al.( 2001a)| compared linear noisy, low-dimensional nonlin- 
ear, stochastic fractional Brownian motion (fBm) and SOC models of magneto- 
spheric fluctuations to examine their scaling and predictabilitity properties. They 
concluded that the SOC model is of particular interest to magnetospheric physics 
due to its robustness in explaining scalings under the wide range of activity levels 
exhibited by the magnetosphere and the solar wind. On the other hand, weakly 
nonlinear models with fBm type noise added could explain the scaling and pre- 
dictability properties of magnetospheric fluctuations equally well. Geomagnetic 
fluctuations on the time scale of substorms and storms certainly appear to be a com- 
pound mixture of multiscale magnetospheric processes including the component 
represented by intermittent solar wind fluctuations [Chang( T999)||Daglis et al.(1999)| . 
Spectral methods making a distinction between these constituents may be queried 
to a certain extent mainly because the information on nonlinear multiscale struc- 
tures is partially hidden for second order statistics. For example, |Tsurutani et al.(1990)1 | 
using the AE-index time series have shown that the power spectrum of the AE data 
exhibits two different power law scalings divided by a spectral break at the fre- 
quency range about l/(5h). While the higher frequency part was thought to be 
more intrinsic to the magnetosphere, the lower frequency part was attributed to the 
influence of solar wind. [ V6rds(2000) | has shown that similar results can be ob- 
tained also by means of a multifractal technique, which also makes possible, how- 
ever, to reveal additional information on singularity distributions of high-latitude 
geomagnetic fluctuations (see later), not evident from spectral studies. Among 
other things, the results of multifractal analysis show that the local singularity 
(Holder) exponents are time dependent and the previously proposed fBm type or 
bicolored noise model [ Tak alo et al.(1993)1 of geomagnetic fluctuations is not rel- 
evant. Moreover, for a range of singularity exponents, a multiplicative cascade 
model, (the P-model) describes quite well the observed singularity distribution 
of geomagnetic fluctuations. Out of this range, however, significant deviations 
from a multiplicative model appear. The P-model describes energy cascade pro- 
cesses in turbulent flows. The largest turbulent eddy is assumed to be built up by 
a specific energy flux per unit length. Then a scale-independent space-averaged 
cascade-rate is considered and the flux density is transferred to the two smaller 
eddies with the same length but different flux probabilities p\ and p^. This pro- 
cess with randomly distributed p\ and p2 (pi + P2 = 1) is repeated towards smaller 
and smaller scales. Energy transfer rate is homogeneous for p = p 1 = p 2 = 0.5 
while p\ > 0.5 corresponds to an intermittent flow. Figure 1 shows how the inter- 
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mittence increases with increasing parameter pi |Tu et al.(1996)| | in a, f(a) plane 
IHals ey et al.(1986")l . The larger spread of a values around right-shifted a, the 
more intermittent field. In spite of the simplicity of the outlined cascade model its 
relevance for the description of intermittence effects in magnetospheric fluctuations 
may lead to the assumption that turbulence rather than F/SOC models fit the ob- 
served statistics better. Some results on statistical distribution of internal time peri- 
ods between bursty geomagnetic events (waiting times) [ K ovacs et al.(20 01)| seem 
to support this assumption. SOC models (at least the original [B ak et al.(1987)| 
model) are expected to display an exponential waiting time distribution. Geo- 
magnetic data, however, display a well-defined power-law waiting time distri- 
bution. It was pointed out by [Boffetta et al.(1999)| that the observed power- 
law waiting time statistics in solar flares appears to be well explained by MHD 
shell models of turbulence. Similar results were presented by |Spada et al.(2001)J 
who analysed density fluctuations in a magnetically confined plasma system and 
found that waiting time statistics is in contrast with the predictions of an SOC 
system. An opposite view was presented by [Free man et al. ( 2000b) | who con- 
jectured that a wider class of running sandpile models | |Hwa and Kardar(1992)| 
could exhibit power law behaviour in the probability density functions of wait- 
ing times. [ Watki ns et al.(2 001b) | determined that the PDFs for burst durations 
and waiting times in a reduced MHD simulation follow power-laws which is not 
sufficient to distinguish between turbulence, SOC-like models and colored noise 
sources. [Boffetta ~al.(1999)| and [Antoni et al.(2001)| also argued that SOC 
models represent self-similar, fractal phenomena. Geomagnetic fluctuations ex- 
hibit clear multifractal scaling [Consolini et al.(1996)| |Voros(2000)| which seems 
to contradict to SOC concepts again. [Georgoulis et al.(1995)|, however, demon- 
strated that the SOC state displayed by their cellular automaton model of isotropic 
and anisotropic energy avalanches has multifractal and multiscaling characteris- 
tics, rather than single power-law scalings and this feature was even enhanced by 
considering extended instability criteria. Similar SOC models of solar flares also 
exhibit multifractal and multiscaling characteristics | Vlahos et al.(1995)| . Recent 
works IVassiliadis et al.(1998)[|lsriker et al.(1998)]|Isliker et al.(2000)||Georgoulis et al.(2001)| 
Urit sky et al.(2001)) underline the proximity between the SOC rules and laws of 
MHD in space physics systems, which are long known to exhibit turbulent be- 
haviour (Georgoulis, personal communication, 2001). More realistic F/SOC mod- 
els which include e.g. less artificial feedback mechanisms, a wide variety of driv- 
ings, interacting avalanches, SOC in continuum physical systems | |Lu (1995)] , etc. 
may further establish the proximity with turbulence. As pointed out by [Urits ky et al.(2001)| 
?] the concept of SOC in a continuum limit (fluid or MHD limit) is essentially un- 
explored. [ Kli mas et al.(200 0) | proposed a simplified Earth's magnetotail current 
sheet model based on continuum SOC model of |Lu (1995)1 - The continuum Lu 
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model in a magnetic field reversal configuration can evolve into SOC due to local- 
ized rapid magnetic field annihilation within the field reversal region. In the same 
time, the plasma sheet is dominated by strong turbulence which keeps the system 
near criticality and produces a predictable quasi-periodic loading-unloading cycle 
of coherent global substorm activity [Klimas et al.(2000)|. In this model turbu- 
lence, SOC states and coherent global modes coexist within the Earth's magnetotail 
on different scales. As a result, the observed ground based and satellite time-series 
contain a mix of fluctuations of different physical origin. [ Angelopoulos et al.(1999) | 
further argued that the presence of intermittent turbulence in the Earth's magneto- 
tail may alter the conductivity and the mass/momentum diffusion properties across 
the plasma sheet and may permit cross-scale coupling processes playing also an 
important role in the establishment of SOC state. 

In this paper no attempt will be made to participate in the theoretical debate on 
SOC or turbulence. However, having a pragmatic view, our opinion is that mea- 
sures of intermittence or char acteristic descriptors of cascade processes, commonly 
used in turbulence studies, like the flatness and multifractal spectra, could be ap- 
plied in solar wind-magnetosphere interaction studies for further comparison of 
basic characteristics of intermittent fluctuations in solar wind and within the mag- 
netosphere. We assume this aproach might be useful providing experimental infor- 
mation on such characteristics of fluctuations which are not accessible for spectral 
studies or second order statistics. This assumption was already investigated in 
IVoros (1998)1 |Voros(2000)]|Voros and Kovacs(2001)l|Kovacs et al(20 01)|. The 
preliminary results allow to make a working hypothesis that intermittence, scaling, 
rapid changes, singularities represent an essential piece of information regarding 
the effectiveness of SW - magnetosphere coupling not considered enough hitherto. 
In order to go deeper, we analyse different geomagnetic and solar wind data sets 
and make a comparison between various magnetic activity levels considering time 
scales of geomagnetic storms (from hours to days), substorms (from half an hour 
to a few hours) or less. 

The main goal of this comparative study is to contribute to the understanding 
of solar wind - magnetosphere interaction processes on the basis of characteris- 
tic scaling/singularity features of the considered time series. We recall that solar 
wind (SW) fluctuations are strongly intermittent [Burlaga(1991) ], that is energy 
at a given scale is not homogeneously distributed in space or/and time. Several 
studies on characteristic probability distribution functions (PDF) of increments 
of SW parameters (magnetic field, velocity, temperature, Elsasser variable, etc.) 
[Mars ch and Tu(1994)||Sorriso-Valvo et al.(19 99) | and on multifractal structure of 



SW fluctuations [Burlaga(1992)[ |Carbone(1994)||Marsch et al.(1996)||Tu et al.(1996)) 



support this assumption. Moreover, [Veltri and Mangeney(1999) | and [Bruno et al.(1999)| 
have shown that there is a direct link between intermittence and the presence of S W 
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structures (Alfvenic, magnetic fluctuations, discontinuities) across which the mag- 
netic field magnitude changes. Also, the high frequency (small scale) fluctuations 
of the southward component of the interplanetary magnetic field have a different 
spectral scaling exponent as the one exhibited by geomagnetic AE index fluctua- 
tions. On larger scales the corresponding spectra are similar [Ts urutani et al.(1 990) | 
A comparative study of dynamical critical scalings in the auroral electrojet (AE) 
index versus solar wind fluctuations confirmed that for times shorter than 3.5 hours 
(higher frequencies) the AE index fluctuations are of internal magnetospheric ori- 
gin |Uritsky et al.(2001)| . 



2 Data analysis methods 

As usual, we introduce the concept of scale (r) through the difference 

6X(t,T)=X(t + r)-X(t) (1) 

where X(t) is the time series under consideration. [Mar sch and Tu(1994)| have 
shown that the PDFs of the increments 6X (X - SW parameters) exhibit strong 
deviations from Gaussianity, especially at smaller scales and the effect is due to 
intermittence of SW fluctuations. To quantify the degree of deviation from the 
Gaussian distribution, i.e. the level of intermittence at different scales, we compute 
the flatness defined by 

F = <^(^) 4 > (2) 

The flatness of a normally distributed signal is equal 3. Adding intermittent fluc- 
tuations to an originally Gaussian signal implies the spreading of its PDF, and 
consequently, the increase of its flatness from the original value of 3. 

Non-homogeneous/intermittent distributions in space/time may also appear as 
asymptotically singular and can be characterized locally, at the point ti, by the sin- 
gularity (Holder) exponents a as [Veh el(1996)| [Vehel an d Vojak(1998)||Riedi(1995)l 
|Canus( 1998)) 

a n {ti) = km (3) 

n— +oo fi 

where \x is a measure constructed from a time series. The procedure involvs the 
computation of the energy content of the differenced signal (Eq. 1) by taking its 
squared value. The measure at a point t{ is given by 5X 2 (ti, r)/ J2i 5X 2 (ti,r). To 
analyse the distribution of singularity exponents (Eq. 3) a sequence of partitions P 
is introduced so that [ Vehel(1996) ] 

P n = {/*}; o < k < 2 n - 1 
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I* =< k2~ n , (k + l)2~ n ) 

where I*; is the interval containing t and the resolution is set by n. The quantity 
of interest is the so-called large deviation singularity spectrum, f(a), which repre- 
sents a rate function measuring the deviation of the observed a from the expected 
value a. The rate function, f(a), can be estimated through [ Vehel(1996) | : 

/(a) = lim (4 ) 

71 — >00 fl 

where N n (a) is the observed number of coarse grain Holder exponents . Usually 
"histogram methods" for estimation of f(a) are used. In that case the number of 
those intervals 1% for which a;„ falls in a box between a m i n and a max is computed 
and /(a) is found by a regression. It yields satisfactory results for pure multi- 
plicative processes, but fails to describe non pure or compound processes when 
f(a) is not a concave function. To overcome this difficulty the so-called double 
kernel method was proposed [Vehel( 1996)1 |Vehel and Vqa k(1998) | realizing that 
N n (a) may be written as a convolution of the density of the a„s and a compactly 
supported kernel. This method allows to estimate non concave rate functions and 
we are going to show that this property may be properly used for characteriza- 
tion of fluctuation processes in near-Earth space. In this paper the estimations of 
/(a) spectra were realized using the FRACLAB package developed at the Institut 
National de Recherche en Informatique, Le Chesnay, France. 

3 Ground based data 

In order to study the basic characteristics of auroral zone geomagnetic fluctuations 
we analyse geomagnetic //-component 1-min mean data from polar cap obser- 
vatory, THULE (THL: 77.47° N, 290.77° E), geomagnetic X-component 1-min 
mean data from high-latitude observatory NARSSARSSUAQ (NAQ: 61.16° N, 
314.57° E) and 1-min mean auroral electrojet (AE) index time series, all from 
1991-1992. 

As known, the AE-index was introduced by [Davis and Sugiura(1966) ] to de- 



scribe the global activity of the auroral zone electric currents and is derived, after 
the substraction of base line values, from evaluation of the variations measured at 
12 stations located near the northern auroral zone. There exist a large number of 
physical mechanisms which couple the auroral zone processes with those within 
magnetospheric tail or in SW. Recently, intermittent energy transport in the mag- 
netotail, the so called bursty bulk flow events (BBF) came into the limelight of the 
magnetosphere research [Baumjohann(1990) , Angel opoulos et al.(1992)| . From 



this point of view the understanding of the response of auroral zone currents or 
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dissipation fields to the time - varying magnetotail dynamics seems to be impor- 
tant. There was not full understanding achieved regarding the nature and origin 
of the related magnetic fluctuations. Partly it was already explained in the intro- 
duction that it is related to the paradigms of F/SOC versus turbulence. We men- 
tion here some other open questions. For example it was shown that the burst 
lifetime distributions of some SW parameters are also of power law form, which 
might be a signature of SOC or turbulence regimes in SW [Freema n et al.(20 00a) , 
Kovacs and Voros (2001) |. Therefore, it is supposed that the scale free property 
of the AE-index may arise from the SW input or at least the internal dynamics 
of the magnetosphere may be masked by the scale free properties of SW driver 
[Free man et al. (2000a) |. Again we remind, however, the very limitations of sec- 
ond order statistics in interpretation of the observed scalings. Another measure of 
the auroral zone dissipation fields is represented by polar optical activity within 
UVI bands. |Lui et al."( 2000) | have examined the blobs of brightness as a proxy 
for BBF events. It was found that the non-substorm "internal" events have a power 
law distribution whereas the system wide events like substorms besides a scale free 
region exhibit a "bump", corresponding to a mean value in substorm breakups. A 
somewhat opposite view was presented by [Conso lini and D e Michelis(1998) | who 
analysed AE-index fluctuations on time scales 1-120 [min] both in quiet (laminar 
phase) and disturbed periods (turbulent phase). They found that in both phases 
the intermittence at different time scales rescales in the same way and the non- 
Gaussian character of PDFs seems to be due to the same physical processes. 

Here we pose the question again about the scaling and singularity proper- 
ties of the AE-index, compared to the similar characteristics of geomagnetic data 
from two observatories THL and NAQ. As the AE index is derived from geo- 
magnetic variations in the horizontal component observed at selected observatories 
along the auroral zone in the northern hemisphere, we expect that the scale r (Eq. 
1) cannot be precisely defined. The geographic distance between observatories 
through Taylor's hypothesis already introduces some effective time shift (scale, r). 
Though, the application of the Taylor hypothesis within the magnetosphere is lim- 
ited [ Dudok de Wit and Krasnoselskikh (1990)| . Besides, from the recordings of 
auroral stations the greatest (upper envelope) and smallest (lower envelope) values 
are taken at intervals of one minute and their difference defines the AE-index. As 
far as the contributing observatory which gives the lower/upper envelope changes 
during the times, and the fluctuations with values between the upper and lower 
envelopes are not taken into account at all (smoothing), the AE-index appears to 
be a measure of auroral zone processes with mixed scales. This is certainly not 
an advantage when a multiscale analysis of AE-index time series is performed. 
This fact is usually neglected in the related literature. A less sophisticated way 
is to take data from a single observatory, but it may have some other drawbacks 
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because significant distant disturbance events can be missed. Also, some BBFs 
with short duration may remain undetected on the ground because of their local- 
ized nature [Daglis ~t al.(1999)| . Nevertheless, a comparison of fluctuations of the 
"multi-observatory measure" (MOM: AE-index) and of the "single observatory 
measure" (SOM: THL geomagnetic field if -component and NAQ geomagnetic 
field X-component) may be instructive. 

Figure 2 shows the PDFs from normalized increments 5X (Eq. 1) of the SOM 
(Figure la, b) and MOM (Figure lb) data (X = X(NAQ), H(THL) and AE, re- 
spectively). Calculations have been made for r = 5, 50, 500, 5000 [min] and 1 
minute mean data was considered from the years 1991-1992. As can be seen the 
distributions in all cases change with r, and for smaller values of r significant de- 
viations from the normal distribution occur. The tails of the distributions reduce 
with increasing scale parameter as a consequence of the decrease of the probabil- 
ity of coherent fluctuations between points separated by increasing distance. The 
difference between the MOM and SOM is more clear if the flatness (F) of the 
corresponding distributions is compared (Eq. 2). Figure 3 shows how the flat- 
ness evolves with increasing r (r £ (5, 5000) [min]). The errorbars correspond to 
the standard deviations (std) computed from time series divided to several parts. 
For small scales, say, r less than a few tens of minutes, the AE-index (MOM) ex- 
hibits smaller deviations from the normal distribution than NAQ (SOM). F(NAQ) 
reaches the level of F max (AE, r = 5) only at the value r ~ 30 [min], which roughly 
may be considered as the effective time shift introduced by the method of deriva- 
tion of the AE-index. THL (polar cap observatory) data essentially show the same 
behaviour as NAQ, but the stds are larger. Therefore, we conjecture that MOMs 
(multi-observatory geomagnetic indices), due to smoothing and scale mixing ef- 
fects, lead to underestimation of the intermittence on small scales. Henceforth, we 
will show only the dependence of the flatness on scale parameter r. 

Let us consider now the singularity distributions f(a) (Eq. 3, 4) for the same 
data sets as above. To see better how the f(a) rate functions evolve with r, all the 
curves are depicted at the same plane (a, /(a)) in Figure 4. Proper symbols are 
introduced for the scales r = 10, 50, 500 [min] (A, *, o) and the dashed lines corre- 
spond to the errors estimated by changing the resolution in Eq. 3, 4. (No averaging 
needs to be done estimating f(a) |Vehel(19"96 ) ] so the spectrum may be evaluated 
at only one resolution. However, to show that the estimations are consistent, all the 
spectra were computed using 15 different resolutions). Again, there are several dif- 
ferencies between the SOM and MOM data. Figure 4a shows the singularity spec- 
tra for NAQ observatory data. On the considered scales, the shape of the curves 
is almost parabolic, close to the P-model fit (thick curve) with p\ ~ 0.745. The 
best correspondence is achieved for the smallest value of r = 10 [min] (depicted by 
symbol A on Figure 4a). It means that in case of auroral zone SOM fluctuations 
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and especially on small scales, the deviations from the Gaussian distribution (Fig- 
ure 3) can be explained by a simple cascade model. Though the phenomenology 
of turbulent cascades in fluid flows is more complex than the simple P-model fit in 
Figure 4a would lead us to indicate, the ID cascade model rouhgly describes how 
the auroral zone SOM fluctuations become more and more intermittent at smaller 
and smaller scales. The spectra for polar cap (THL observatory) SOM (Figure 4b) 
and AE index MOM (Figure 4c) have a more pronounced non-parabolic shape indi- 
cating the presence of compound processes. At small scales (r = 5,50 [min]) THL 
observatory fluctuations contain stronger singularities because of the extension of 
the rate functions left wings to the smaller (more singular) values of a. For r = 
500 [min], however, the right wing evolves to the less singular values which may 
be related to the SW influence [Voros(2000)|. This effect is less visible, but still 
present in auroral zone SOM data (Figure 4a). The MOM AE index /(a) spectra 
practically do not change with r. We conjecture, this is the result of the method of 
derivation of the AE index, resulting smoothing and scale mixing effects. 

In all cases, the deviations of singularity spectra from the parabolic shape may 
be indicative of the phenomenon of phase transition. Namely, at the a values where 
the /(a) spectra are out of parabolic shape the major contributor to the observed 
singularities may change from one measure to another. As different physical pro- 
cesses may generate different measures (distributions), possible models with simi- 
lar characteristics as the observed spectra may contain physical information on the 
contributing (e.g. SW or magnetospheric) sources. 

4 Satellite data 

The very advantage of the ground based data is its availability for a long periods of 
time. For a proper estimation of PDFs or singularity rate functions long data sets 
are needed which is a requirement hardly ever fulfilled in the case of satellite data. 
Nevertheless, we expect to find out some interesting scaling/singularity features of 
interplanetary magnetic field (IMF) fluctuations proceeding in the same way as in 
the previous section. To this end we analyse ACE and WIND IMF B magnitude 
and B z component data which are available with time resolution of 16 [s] and 
3 [s], respectively. SW velocity is not considered here due to too many gaps in 
data. While the ACE satellite is continuously monitoring the SW at the L\ point, 
WIND has a more complicated trajectory crossing also the magnetosphere from 
time to time. For our analysis we have chosen time periods when WIND was also 
in SW and there were negligible data gaps in both cases (less than 1 % of total data 
lengths). 

Another aim was to analyse "geoeffectively different" periods of IMF B and 
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B z magnitude fluctuations. Geoeffectiveness during the chosen periods was con- 
sidered examining the geomagnetic D st index which is derived from the geomag- 
netic field //-component registrations of 4 observatories [Sugiura(1964) | and it 
aims at giving the effect of the magnetospheric ring currents. The chosen peri- 
ods were classified as disturbed ones if geomagnetic fluctuations with storm-index 
D st being less than -50 [nT] occured several times within a consedered interval. 
The limit of D st < -100 [nT] which corresponds to intense storms was con- 
sidered, too. We emphasize, however, instead of a study of individual storms, 
the generic features of fluctuations on a given scale r, but during longer periods 
of time are investigated. Essentially 2-^4 weeks of data with the consedered 
time resolutions may already ensure sufficiently robust estimations of singular- 
ity spectra (Eq. 4). In this sense, several intense magnetic storms D st -C -50 
[nT] may occur during a strongly disturbed interval, a less disturbed period con- 
tains less intense storms and an undisturbed period has only D st > -50 [nT]. The 
limit of -50 [nT] was chosen on the basis of previous studies of magnetic storms 
[Tayl or et al.(1996)| . Intense magnetic storms are characterised by D st index < 
-100 [nT] |Gonzalez and Tsurutani(1987)| . In this preliminary study of generic 
features of magnetic fluctuations these limits may be consedered as more or less 
adequate. We expect that this rough classification of the geomagnetic response al- 
lows us to identify characteristic scaling and singularity features of the correspond- 
ing IMF magnetic fluctuations that would be indicative for their geoeffectiveness. 

There were 5 time periods and 6 data sets separated for our analysis (for one 
period there were both ACE and WIND data available). The disturbed periods 
are the following: 1. March 19 - April 25, 2001 (ACE); 2. October 1 - November 
30, 2000 (ACE); 3. April 9 - April 20, 1997 (WIND). The undisturbed periods: 4. 
November 18 - December 10, 1998 (ACE); 5. January 10 - January 29, 1998 (ACE, 
WIND). For demonstration we show some of the data sets. 

Figure 5 shows the first ACE data set from March 19 to April 25, 2001. For 
the sake of perspicuity, the time resolution is 1 [h]. (The flatness and singularity 
spectra are computed from the time series with time resolution of 16 [s] and 3 [s].) 
During this extremely active period intense magnetic storms (D st < -100 [nT]) 
occured several times (Figure 5c). The limit of -100 [nT] is depicted by a thick line 
in Figure 5c. [Gonzale z and Tsurutani(1 987) [ have shown that the interplanetary 
causes of intense magnetic storms are long duration (> 3 [h]), large and negative 
(< -10 [nT]) IMF B z events associated with interplanetary duskward electric fields 
> 5 [mVm^ 1 ]. In Figure 5b IMF B z is depicted, including a thick line indicating 
the level of -10 [nT]. Comparison of Figures 5b, c shows an agreement with the 
above criteria, that is, long duration negative IMF B z events occur together with 
intense magnetic storms. Figure 5a shows the variations of IMF B. It is visible 
that an intense magnetic storm occured at the end of the studied period, between 
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t = 800 and 900 [h] (Figure 5c). B z < 10 [nT] (Figure 5b) and B max ~ 17 [nT] 
corresponded to this event. A similar enhancement of B at 400 [h] < t < 450 
[h] appears in Figure 5a having no intense storm response in D st , which can be 
explained by the corresponding IMF B z > -10 [nT] in Figure 5b. 

Figure 6 shows an undisturbed period from January 10 to January 29, 1998. 
Again, IMF B, B z (WIND, ACE Figures 6a-d) and D st index (Figure 6e) are 
shown. It is visible that D st > -50 [nT] and B z > -10 [nT] everywhere. 

Fluctuations of IMF and their geoeffectiveness were studied by a number of 
authors. [McPherr on et al.(1986)| showed that substorms are frequently triggered 
by changes in the IMF. [Kamide(2001) | proposed that the quasi-steady compo- 
nent of the interplanetary electric field is imporant in enhancing the ring current, 
while its fluctuations are responsible for initiating magneospheric substorms. It 
is out of scope of this paper to analyse the influence of other interplanetary pa- 
rameters (e.g. velocity, density, temperature, etc.) on storm/ substorm activity 
[Dagl is et al.(2001)| . Rather we will concentrate on the level of intermittence of 
IMF fluctuations. 

To this end let us consider the flatnesses and the singularity rate functions for 
the disturbed and undisturbed periods of 1 — 5. Deviations from the Gaussian 
distribution are larger in the case of disturbed events depicted by larger markersizes 
in Figure 7. The PDFs of the undisturbed events are also non-Gaussian, but the 
flatnesses for a given scale are smaller than those for disturbed events especially at 
scales t < 1000 [s]. Similar differences are present in singularity spectra computed 
for IMF B fluctuations at the scales of r = 320 [sec] (ACE) and r = 60 [s] (WIND) 
shown in Figure 8. There were the same marker types used as in Figure 7. ACE and 
WIND data are depicted separately in Figures 8a, b. The maxima of the singularity 
spectra of more disturbed periods have a tendency for shifting to larger values of 
a. Also, the spread of singularities around most probable a is wider for more 
disturbed cases. But it is the same behaviour as in case of the simple P-model in 
Figure 1, when intermittence is stronger and stronger for larger and larger values of 
pi. The differences between the disturbed and undisturbed cases gradually cease 
for larger values of r (not shown). 

Figure 9 shows the singularity spectra computed for IMF B z fluctuations in 
the same way as previously. Obviously, the intermittent fluctuations of IMF B and 
B z fields exhibit very similar changes in their flatnesses and singularity spectra as 
the geoeffectivity level changes. For example, there is a clear difference between 
the introduced scaling and singularity characteristics (Figure 5) for the disturbed 
period March 19 - April 25, 2001 and for the undisturbed period January 10 - 
January 29, 1998 (Figure 6). It indicates that, in addition to known "geoeffective" 
SW parameters (e.g. southward component of IMF) or their combinations, small 
scale rapid changes, singularities and non-Gaussian statistics of IMF fluctuations 
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may play an important role in SW - magnetosphere interaction processes. 

Also, the question naturally arises to what extent the magnetospheric response 
itself is influenced by small scale statistics of IMF fluctuations. Unfortunately, not 
all the geomagnetic data is available for the above analysed periods. It is possi- 
ble to test, however, how the shape of the rate function changes if available data 
is considered. Minute-mean i7-component geomagnetic data from THL observa- 
tory is available from 1975 to 1996. We computed the f(a) spectra for each year 
at the scale r = 50 [min] and analysed how their shapes change at the f(a) val- 
ues 0.3 - 0.9. At each level of f(a) the corresponding values of a m i n and a m ax 
were computed (filled circles in Figure 10a). In Figure 10b the time evolution 
of the difference a max — a m i n at a given f(a) level is depicted. The average 
standard deviation at f(a) = 0.3 is about 0.06, while at f(a) = 0.9 is 0.02. At 
f(a) = 0.3 - 0.6 the a max — a m i n curves strongly fluctuate indicating significant 
changes in the shape of f(a) spectra from 1975 to 1996. Figure 11 shows the 
similar results of [Kam ide et al.(1998 ) |, however, obtained by a different method. 
[ Kam ide et al.(1998 ) ] analysed the occurence of geomagnetic storms in compari- 
son with yearly averaged Wolf sunspot number. Figure 1 1 shows the yearly aver- 
aged number of hours with D st less than -100 [nT] (solid line with filled diamonds), 
and with D st less than -50 [nT] (divided by 5, dashed line). Thick line corresponds 
to yearly averaged sunspot number. It is visible that the maxima of geomagnetic 
activity and of solar cycle do not coincide. During the declining phase of the so- 
lar cycle coronal holes emerge from polar regions of the Sun which are continuous 
sources of fast-speed plasma causing a peak in recurrent geomagnetic storms activ- 
ity [Kamide et al.(1998) Kamide(2001) |. The similarity between the variability of 
rate function shapes for f(a) = 0.3 0.6 (Figure 10b) and yearly averaged number 
of hours with prescribed D st indices (Figure 11) is remarkable. This correspon- 
dence also supports our working assumption that the shape of /(a) rate function 
estimated using Eqs. 3, 4 [Veh el(1996)||Vehel and Vqj ak(1998) | contains relevant 
physical information. We mention that available AE-index and NAQ observatory 
data lead to the same results. On the other hand, however, at /(a) = 0.7 - 0.9 the 
variations of a max — a m i n versus time are negligible (Figure 10b). 

On the basis of Figures 10b and 1 1, years with extreme levels of geomagnetic 
response can be chosen and their corresponding singularity spectra can be recalcu- 
lated. Figures 12a-c show the SOM and MOM spectra estimated at r = 50 [min] 
for maximum (1991, symbol o) and minimum (1984 - symbol x) years of geomag- 
netic activity. One can see that the characteristic asymmetric shape of AE-index 
spectra (Figure 12c) in comparison with Figure 4c is present henceforward, pre- 
sumably caused by the influence of SW fluctuations. Noticeably, the maxima and 
the right wings of the more geoeffective SW IMF singularity distributions (Figures 
8, 9) match well the region of a S (1.2, 1.8) within which the MOM spectra for 
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the most part are out of parabolic shapes. The same effect is visible for the more 
singular wing of f(a) spectra at a G (0.4, 0.8). The asymmetry is even enhanced 
in 1991 (maximum of geomagnetic activity). The differences between more active 
(1991) and less active (1984) years are present in SOM spectra mainly at the wings 
of /(a) rate function, too (Figures 12a, b). The most probable singularities around 
a do not exhibit any changes. 

5 Discussion 

We presented an analysis of scaling and singularity characteristics of ground based 
and satellite magnetic fluctuations in this paper. These techniques are commonly 
used in studies of turbulent flows. As we oulined in the Introduction the sup- 
posed proximity between F/SOC and MHD turbulence models may allow to esti- 
mate measures of scaling, intermittence, etc., directly from time series, for futrher 
comparison aiming to help a development of tractable numerical models of highly 
variable SW - magnetosphere interaction (Georgoulis, personal communication, 
2001). We think the results are not contradictory in this sense, but rather eluci- 
date important aspects of magnetic field fluctuations which should be incorporated 
to more realistic F/SOC models of SW driven magnetospheric activity. Scaling 
and singularity characteristics of high-quality geomagnetic data from THULE and 
NARSSARSSUAQ observatories (Single Observatory Measures) and of AE-index 
(Multi Observatory Measure) were compared estimating PDFs, flatnesses and sin- 
gularity rate functions. The same methods were applied for ACE and WIND data. 
The specific periods chosen for the analysis of SW fluctuations reflect the limited 
availability of high resolution satellite data for relatively longer periods of time. In 
spite of this, subgroups of disturbed and undisturbed periods were selected, having 
primarily in mind the occurence of enhanced geomagnetic response represented 
by 1-hour mean D st -index. The comparison of generic features of fluctuations in 
SW and high-latitude SOM and MOM data is then possible, because geomagnetic 
activity at high-latitudes is always very high during magnetic storms, though the 
storm/substorm relationship itself is more complicated [Dagli s et al.(20 01)|. 

It was shown that the intermittence of AE-index fluctuations is reduced at small 
scales due to the method of its derivation. We argue this also provides a possible 
explanation for the negligible changes of the AE-index singularity distribution with 
r. Other kind of MOM data might be influenced in the same way. Nevertheless, 
the contribution of SW fluctuations make the AE-index rate function asymmetric, 
mainly within the range of a G (1.2,1.8). The same asymmetry caused by the SW 
is also present in THL and NAQ rate functions, however, mainly for r > 300 [min] 
(see also |Voros(2000)l ). 
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In case of SW fluctuations it was demonstrated that the departure from Gaus- 
sianity is stronger at small scales and it clearly depends on the correlation (geo- 
effectiveness) between IMF fluctuations and the occurence of geomagnetic storms 
(decreased D st -index). It seems to indicate that the intermittence strength of IMF 
magnitude and B z component fluctuations, in addition to other SW parameters 
such as southward B z , SW velocity, density, Alfvenic Mach number, plasma j3, 
represents a new parameter (or rather a whole set of parameters describing singu- 
larity features) controlling the energy input rate to the magnetosphere. Considering 
Taylor's hypothesis and SW velocities of 500 km/s, intermittence at time scales of 
tens of seconds corresponds to the spatial structures of several thousand kilome- 
ters or more. |Book and Sibeck( 1995) 1 estimated the corresponding timescale on 
which turbulent motion may affect the transport of mass and energy across the 
magnetopause through interchange instability. They have found it is less than 150 
[s]. As known from previous ISEE1 and 2 magnetometer studies day side recon- 
nection of IMF and GMF lines seems often to be a sporadic and patchy process 
and measurements obtained at and near the magnetopause indicate that reconnec- 
tion does not necessarily occur across the all dayside magnetopause even under 
the favourable southward pointing IMF conditions [Rijnbee k et al.(19 84) [. We 
conjecture, patchy reconnection may be related to intermittence and singularity 
characteristics of IMF turbulence at small scales. There is a number of works in 
which the role of turbulence in magnetopause reconnection processes is anticipated 
IGaleev et al.(1986)| |Drake et al.(1994)l |Kuznetsova and Roth(1995)| . We recall 
the work of [ Galee v~et al.(1986)"| in which patchy reconnection was considered to 
be an irregular multiscale process associated with the magnetic field diffusion and 
self-consistently generated magnetic turbulence. Our results indicate that a num- 
ber of singularity parameters (Holder exponents) should be taken into account to 
properly describe the basic characteristics of the upstream SW turbulence. In this 
paper we examined the global distribution of IMF singularities and found clear 
differences between geoeffectively disturbed and undisturbed periods. Obviously, 
to understand better the role of turbulence in patchy magnetopause reconnection 
processes a proper time and space localization of IMF singularities will be needed. 

As far as the magnetospheric response is considered, previous results |V oros and K ovacs( 2001)[ 
have suggested that global singularity spectra estimations of SOM and MOM data 
sets on different scales may allow to separate fluctuations of SW or magnetospheric 
origin. Our results show that the influence of the SW is perceptible mainly at the 
wings of the rate function, that is at smaller values of f(a). The most probable 
singularities (f(a) = 0.8 — 1) are less influenced by SW driver. Rate functions 
estimated for years 1975-1996 exhibit similar variations as geomagnetic activity 
studied by [Kamide et al.(1998)|. SW forcing effects were found when SOM and 
MOM singularity spectra for two years (1984 and 1991) of different geomagnetic 
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activity levels were compared. 

We believe that further development in this direction will result a better under- 
standing of SW - magnetosphere interaction allowing more efficient prediction of 
space weather. 
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The P model 
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Figure 1 : The P-model for three different values of parameter p\ 
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Figure 2: Probability distribution functions computed for: a. NARSSARSSUAQ 
X-component geomagnetic field; b. THULE //-component geomagnetic field; c. 
AE-index data. The time scales are r = 5, 50, 500, 5000 [min]. Continuous lines 
refer to Gaussian distribution, while the "o" markers correspond to the analysed 
data. 
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Figure 3: The flatness as a function of scale (r G (5, 5000) [min]) estimated for 
NARSSARSSUAQ (A), THULE data (o) and AE-index (•). 
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Figure 4: Dependence of singularity rate functions on scale; r = 
5(A);50(*);500(o) [min]. a. NARSSARSSUAQ data, b. THULE data; c. AE- 
index. Thick line in Figure 4a corresponds to P-model fit. 
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Figure 5: An extremely disturbed period as seen in 1-hour mean a. ACE IMF 
magnitude, b. ACE IMF B z component and c. D s t-'mAex data plots. 
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Figure 6: Undisturbed period as seen in 1-hour mean a. IMF B (WIND); b. IMF 
B z (WIND); c. IMF B (ACE); d. IMF B z (ACE) data plots and e. the correspond- 
ing geomagnetic response represented by D st -mdex. 
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Figure 7: The flatness as a function of scale (r <G (16*5,16*10 5 ) [sec]) estimated 
for the periods: March 19 - April 25, 2001, (ACE: A); October 01 - November 30, 
2000, (ACE: o); April 09 - April 20, 1997, (WIND: *); November 18 - December 
10, 1998, (ACE: y); January 10 - January 29, 1998, (WIND: x); January 10 - 
January 29, 1998, (ACE: o). Disturbed, that is geoeffective events are depicted by 
larger markersizes. a. IMF B, b. IMF B z . 
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Figure 8: Singularity spectra (rate functions) estimated for IMF B. The periods 
and data sets are depicted by the same marker types as in Figure 7. 
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Figure 9: Singularity spectra estimated for IMF B z . The periods and data sets are 
depicted by the same marker types as in Figure 7. 
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Figure 10: Time variation of the shape of singularity spectra a. filled circles show 
the points in which a m i n and a max values are computed; b. a max — a m i n values 
computed at different f(a) levels. 




Years 

Figure 11: Yearly averaged number of hours with D st < -100 nT (o), and with D st 
< -50 nT (— — — ). Yearly averaged sunspot number (thick line) is also shown. 
After |Kamide et al.(1998)| . 
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YEARS: 1984 (x); 1991 (o) 
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Figure 12: Singularity spectra for a. NARSSARSSUAQ data; b. THULE data; c. 
AE-index. Maximum of geomagnetic activity in 1991 (o); minimum of geomag- 
netic activity in 1984 (x). 
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